DUDLEY KHOX LIB" ' RY 

naval pc ' .l:sc;l U l 
MOMT' iicY. c.._: rora.'iA 93943 



NAVAL POSTGRADUATE SCHOOL 

Monterey, California 




THESIS 

THE EFFECT OF INTERIOR MOTION 
ON 

SEASONAL THERMOCLINE EVOLUTION 
by 

Janice P. Garner 
December 1983 

Thesis Advisor: R. W. Garwood, Jr. 

Approved for public release; distribution unlimited 






Unclassified 



SECURITY CLASSIFICATION of THIS PAGE <TW>»n D«« Enter ea) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT’S CATALOG NUMBER 


4. TITLE (and Subtllla) 

The Effect of Interior Motion on 
Seasonal Thermocline Evolution 


5. TYPE OF REPORT & PERIOD COVERED 

Master's Thesis 
December 1983 


6. PERFORMING ORG. REPORT NUMBER 


7. AUTHORfF) 

Janice P. Garner 


8. CONTRACT OR GRANT NUMBERftJ 


*■ performing ORGANIZATION name ano address 

Naval Postgraduate School 
Monterey, California 93943 


10. PROGRAM ELEMENT. PROJECT, TASK 
AREA A WORK UNIT NUMBERS 


tl. CONTROLLING OFFICE NAME ANO AOORESS 

Naval Postgraduate School 
Monterey, California 93943 


12. REPORT OATE 

December 1983 


13. NUMBER OF PAGES 

85 


14. MONITORING AGENCY name A AOORESS (It different from Controlling Office) 


IS. SECURITY CLASS, (o 1 thla report) 

Unc las s if ied 


15«. DECLASSI FI CATION/ DOWN GRADING 
SCHEDULE 



l«. DISTRIBUTION STATEMENT (ot thte Rapott) 

Approved for public release, distribution unlimited. 



17. DISTRIBUTION STATEMENT (of the ebetrect entered In Block 20, It different from Report) 



IS. SUPPLEMENTARY NOTES 



IS. KEY WORDS (Continue on reveree mldo If neceeeery and Identity by block number) 

mixed layer modelling, vertical advection, 
long period interior motion. 

SO. ABSTRACT (Continue on reveree aide It neceeeery end Identify by block number) 

The response of the seasonal thermocline formation and 
mixing to prescribed vertical interior motion is examined. 
For the annual and shorter period interior motion cases 
studied, the response was strongest for the annual period. 
For an oscillatory vertical motion having a 15 m amplitude 

Unclassified 



DO | JAN*7S 1473 EDITION OF I NOV «S IS OBSOLETE 
S/N 0102- LF- 014- 6601 



SECURITY CLASSIFICATION OF THIS PAGE fl*7>.n Data gniara , 



Unclassified 

SECURITY CLASS! TIC ATI ON OF THIS PAGE (Whn Dmtm Enfrmd) 



at a 75 m depth, warm sea surface temperature anomalies of 
up to 2.16 C and cold anomalies of up to 2.04 C are found, 
depending on the phase difference between the interior motion 
and the annual heat cycle. The net phase-averaged effect 
on the mixed layer for annual period interior motion is a 
reduction in vertical mixing. Higher frequency motion 
produces a net enhancement of the mixing. 



$ N 0102- LF. 014- 6601 



Unclassified 



SECURITY CLASSIFICATION OF THIS P AGEfWh»n Dmtm Entmrmd) 



Approved for public release; distribution unlimits 



The Effect of Interior Hotion 
on 

Seasonal Thermocline Evolution 



by 



Janice P. ^Garner 
Lieutenant, United States Navy 
B. A « , Montclair State College, 1972 
M.A., Indiana University, 1974 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN METEOROLOGY AND OCEANOGRAPHY 



from the 

NAVAL POSTGRADUATE SCHOOL 
December 1983 



f 




DUDLEY KrGY LI r '~ 1 
NAV’L ? r ' 
y- ■ /. 




ABSTRACT 



The response of the seasonal thermocline formation and 
mixing tc prescribed vertical interior motion is examined. 
For the annual and shorter period interior motion cases 
studied, the response was strongest for the annual period. 
For an oscillatory vertical motion having a 15 m amplitude 
at a 75 m depth, warm sea surface temperature anomalies of 
up to 2.16°C and cold anomalies of up to 2. 04 °C are found, 
depending on the phase difference between the interior 
motion and the annual heat cycle. The net phase-averaged 
effect on the mixed layer for annual period interior motion 
is a reduction in vertical mixing. Higher frequency motion 
produces a net enhancement of the mixing. 



4 



TABLE OF CONTENTS 



I. INTRODUCTION 11 

II. MODEL AND EQUATIONS 15 

A. THE MODEL 15 

1. Seme General Considerations 15 

2. Equation of State 15 

3. THE HEAT EQUATION 16 

1. Initial Form 16 

2. Scale Analysis and Simplification of 

the Heat Equation 20 

C. MIXED LAYER DEPTH 21 

D. INTERIOR VERTICAL MOTION 23 

E. POTENTIAL ENERGY AND HEAT 27 

III. EXPERIMENTS AND RESULTS 32 

A. SUMMARY OF PERTINENT EQUATIONS 32 

3. RESULTS WITH NO INTERIOR MOTION 34 

C. ONE-YEAR PERIOD INTERIOR MOTION 36 

1. Effect on Mixed Layer Depth 37 

2. Effect on Mixed Layer Temperature .... 39 

3. Effect on Heat Content 41 

4. Effect on Potential Energy 44 

5. Summary of Response to One-Year Period 

Interior Motion 45 

D. EFFECT OF SHORTER PERIOD INTERIOR MOTION ... 47 

E. ASYMMETRY IN ENHANCED AND REDUCED MIXING ... 48 

IV. CONCLUSIONS AND RECOMMENDATIONS 52 

A. CONCLUSIONS 52 

B. RECOMMENDATIONS 52 

5 



APPENDIX A: FIGURES 55 

BIBLIOGRAPHY 82 

INITIAL DISTRIBUTION LIST 83 



6 



LIST OF FIGURES 



2.1 Turbulent Heat Flax Profile 13 

2.2 Initial Temperature Profile 19 

3.1 Net Surface Heat Flux 33 

3.2 Heat Content Example 1 42 

3.3 Heat Content Example 2 43 

3.4 Two-Layer Ocean Example 50 

A. 1 Temperature, No Interior Motion 55 

A. 2 Temperature, One-Year Period, 90° Phase .... 56 

A. 3 Temperature, One-Year Period, 270° Phase .... 56 

A. 4 Mixed Layer Depth, Nc Interior Motion 57 

A. 5 Mixed Layer Temperature, No Interior Motion . . 57 

A. 6 Heat, No Interior Motion 58 

A. 7 Potential Energy, No Interior Motion 58 

A. 8 One-Year Period Interior Motion 59 

A. 9 h (meters) , One-Year Period Interior Motion . . 60 

A. 10 Ah (meters). One- Year Period Interior 

Morion 60 



A. 1 1 T„ l (C) , One-Year Period Interior Motion .... 61 

A. 12 AT Ml (C) , One-Year Period Interior Motion ... 61 

A. 13 T MU (C) , One-Year Period, Second Year 62 

A. 14 AT mu (C) , Cne-Year Period, Second Year .... 62 

A. 15 H (100 *C m) , One-Year Period Interior Motion . . 63 



A. 16 AH (100 ®c m) , One-Year Period Interior 

Motion 63 

A. 17 p (10 4 °C m 2 ), One-Year Period Interior 

Motion 64 

A. 18 AP ( 1 0 4 °C m 2 ) , One-Year Period Interior 

Motion 64 



7 



A. 19 1/2-Year Period Interior Potior. 65 

A. 20 h (meters) , 1/2-Year Period Interior Motion . . 66 

A. 21 Ah (meters), 1/2- Year Period Interior 

Motion 66 

A. 22 T M1 _ (C) t 1/2-Year Period Interior Motion .... 67 

A. 23 (C) , 1/2-Year Period Interior Motion ... 67 

A. 24 T ML (C) , 1/2-Year Period, Second Year 68 

A. 25 AT ml (C) , 1/2-Year Period, Second Year .... 68 



A. 26 H (100 °C m) , 1/2-Year Period Interior Motion . . 69 



A. 27 AH (100°C m) , 1/2-Year Period Interior 

Motion 6 9 

A. 28 P (10 4 ®C m 2 ), 1/2- Year Period Interior 

Motion 70 

A. 29 AP ( 1 0 4 m 2 ) , 1/2-Year Period Interior 

Motion 7 0 

A. 30 1/3-Year Period Interior Motion 71 

A. 31 h (meters), 1/3-Year Period Interior Motion . . 72 

A. 32 Ah (meters), 1/3- Year Period Interior 

Motion 72 

A. 33 T HL (C) , 1/3-Year Period Interior Motion .... 73 



A. 34 AT ML (C) , 1/3-Ysar Period Interior Motion ... 73 

A. 35 H (100 °C m) , 1/3-Year Period Interior Motion . . 74 



A. 36 AH (100 a C m) , 1/3-Year Period Interior 

Motion 74 

A. 37 P ( 10 4 ° C m 2 ), 1/3- Year Period Interior 

Motion 75 

A. 38 AP (10 4 °C m 2 ), 1/3-Year Period Interior 

Motion 75 

A. 39 1/4-Year Period Interior Motion 76 

A. 40 h (meters), 1/4-Year Period Interior Motion . . 77 

A. 41 Ah (meters), 1/4-Year Period Interior 

Motion 77 



A. 42 T ML (C), 1/4-Year Period Interior Motion .... 78 



8 



A. 43 A T «-il ( c ) • 1/4-Year Pariod Innericr Motion ... 78 

A. 44 H (100 °C m) , 1/4-Year Pariod Interior Motion . . 79 



A. 45 AH (100 C m) , 1/4-Year Period Interior 

Motion 79 

A. 46 P (10* °C m 2 ) , 1/4- Year Period Interior 

Motion 80 

A. 47 AF (10 4O C m 2 ) , 1/4-Year Period Interior 

Motion 80 

A. 48 Peak Enhanced and Reduced Cooling Values, 

Year 1 81 

A. 49 Peak Enhanced and Reduced Cooling Values, 

Year 2 81 



9 



ACKNOWLEDGHENT 



I wish to thank Professors R. W. Garwood and A.J. 
Willmctt for their valuable help ana encouragement during 
the undertaking of this thesis. 



10 



I. INTRODUCTION 



Interior vertical motion can play a significant role in 
the mechanical energy budget of the ocean boundary layer. 
Linden (1975) carried out experiments to observe mixed layer 
deepening in a stratified fluid. He found that internal 
waves are important in the energetics of the mixed layer and 
that the mixing process may provide a significant source of 
energy for the generation of infernal waves. In this thesis 
the interior motion is prescribed and not affected by the 
mixed layer processes. However, the modelling of the inter- 
action of free internal waves with mixed layer processes is 
an important area for future study. Linden points out that 
he did not account for waves generated elsewhere which 
propagate into the region of consideration. Also, the 
internal waves considered by Linden were of relatively high 
frequency. 

De Szoeke (1960) found wind stress curl-driven vertical 
advecticn to be as important as the entrainment rate in 
mixed layer dynamics. Cus hman-Roisin (1981) likewise found 
wind-driven advecticn to be as important as wind-driven 
mixing in mixed layer dynamics and temperature front 
formation. 

Much work has been done to study the ocean response to 
storms. Price et al. (1978) , in their work or. the mixed 
layer response to stcrms, found it necessary to include 
vertical advection. Price (1981) found that upweiling will 
significantly enhance the sea-surface temperature response 
to a slowly moving hurricane. Adamac at al. (196 1), using a 
two-dimensional model to study the oceanic response to a 
hurricane, found a strong interdependence between mixing and 
advection. Sreatbatch (1983) also found advecticn 



to be 



important in increasing maximum cooling due no the passage 
of a storm. 

Only recently has much been done to include long period 
interior motion in mixed layer models. Stevenson 

(1980,1981) examined the response of the mixed layer when 
perturbed by linear Eossby waves. He found no mixed layer 
response during the first half of the heating season while 
the mixed layer is shallowing and a maximum sea surface 
temperature response at the end of the heating season and 
beginning of the cooling season. Stevenson found that the 
maximum response of mixed layer depth to the wave occurred 
at or near the end of the cooling season. Also, he found 
that the lowest frequency waves produced the greatest mixed 
layer depth and temperature responses. This study differs 
from that of Stevenson in that a finite amplitude prescribed 
interior motion is used rather than a linear Rossby wave. 
This produces a phase dependent asymmetry in the mixed layer 
response to the interior motion. 

Burger (1 982) . successfully introduced tidal period 
vertical advection into the Garwood (1976,1 977) one- 
dimensional mixed layer model. He found that the interior 
motion had a significant effect on the near-s urf ace tempera- 
ture field and that the inclusion of the interior motion can 
improve the model for single station forecasting. 

In this study, the Garwood (1976,1977) one-dimensional 
ocean planetary boundary layer model is used. Vertical 

advection was added to the model. The annual heat cycle is 
a sine function and wind stress is held constant at one dyne 
per square centimeter. By keeping the surface boundary 

conditions simple, the effect of the interior motion can be 
more readily observed. 

The interior vertical motion is a prescribed, forced 
interior wave which is linear with depth and sinusoidal in 
time : 
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( 1 . 1 ) 



W ( 2. , t ) 
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where D is a reference depth, 200 meters, and W 0 is the 
amplitude cf the prescribed vertical motion at depth, D. 
The periods, T, studied are one year, one-half year, one- 
third year and one-quarter year. The phase difference, <f> , 
between the interior motion and the annual heating cycle is 
varied from 0° to 360° at increments of 22.5°. This 
interior motion is general and may be representative of a 
number cf possible long period phenomena, for example, 
internal Ecssby waves or perhaps Ekman pumping on seasonal 
scales. 

The oceanic variables of interest are mixed layer depth, 
h, sea surface temperature which equals mixed layer tempera- 
ture, T mi _, the heat content, H, and the potential energy, P. 
These variables are contoured as functions of time and phase 
for each cf the pericds studied. Also, contours are made of 
the differences of these variables relative to the case of 
no interior motion. ‘These contours are used to observe the 
effect of the interior motion on the mixed layer processes. 

Furr her mere, the phase, <? , in equation (1.1) can be 
replaced by the expression: 

4 = , ( 1 . 2 ) 

L* 



where x is horizontal dis placemen- and L is some typical 
wavelength. Then, at any fixed x, equation (1.2) represents 
the vertical motion due to the passage cf a travelling wave. 
The wavelength, L, might be, for example, 300 km for an 
annual karoclinic Eossby wave in the central north Pacific 
(Kang and Kagaard, 1980). The phase in equation (1.1) can 
thus be related to horizontal location. The contours of h, 
T mi _ , H and P versus phase and time then indicate horizontal 
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structure in these 



variables due 



to the presence 



wave . 



of the 
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II. MODEL AND EQUATIONS 



A. THE MODEL 

1 . Some Ge ner al Con sid erati ons 

For this study the Garwood (1976,1977) one- 
dimensional bulk mixed layer model is used. The model is 
based on the Na vier-Stokes equations of motion, assumes the 
continuity equation for an incompressible fluid and makes 
the one-dimensional approximation. The Garwccd model 

differs from previous bulk mixed layer models in several 
features. The time scale for viscous 'dissipation cf turbu- 
lent kinetic energy, TKE, is dependent on planetary rota- 
tion. This allows for a cyclical steady state on an annual 
basis by increasing dissipation for deeper mixed layers. 
Eoth downward and upward turbulent heat fluxes are important 
in this model. Furthermore, entrainment depends on the 
relative distribution of TKE between vertical and horizontal 
components. 

Before continuing, a comment on notation is in 
order. An overbar cn a variable will indicate the horizon- 
tally averaged mean. A prime on a variable will indicate 
the fluctuation about this horizontal mean. Angled 
brackets, < >, will indicate vertically averaged, mixed 

layer mean. 

2. Equ ati on of S tat e 

In the Garwood model, the equation of state has the 

form : 
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( 2 . 1 ) 



r = fo [ I - * (T-T 0 ) - p(S- 5- 0 )] , 

The thermal expansion coefficient, ^ , is a function of T 

and is re-computed daily. The constants .a., T 0 and are 
typical density, temperature and salinity values, respec- 
tively. In this study, salinity is fixed (precipitation and 
evaporation are disallowed) and set equal to the value of 
S p . Equation (2.1) then becomes 

p = /’o[l - ^(T-T„)3 . < 2 - 2 ) 



B. THE HEAT EQUATION 
1 • Ini tia l Form 



The time rate of change of temperature is given by 
the one-dimensional heat equation: 



^ i 
o>t 




JlL 




w IX 

& 



-+ K 






(2.3) 



The first term on the right is vertical turbulent heat flux, 
the second is solar insolation, the third is vertical 
temperature advection and the last, diffusion. All of the 
solar insolation will be assumed to be absorbed at the 
surface (in the first few centimeters) of the ocean. Hence, 
the solar insolation term will be dropped in favor of a 
surface boundary condition on the turbulent heat flux. 
Equation (2.3) then becomes 




(2.4) 
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mi xed 



In this model , 
to be vertically homogeneous, 
to 



layer temperature is assumed 
Equation (2.4) then reduces 



Tf.lL 

}-k 




w'T‘ ) 



(2.5) 



for mixed layer temperature, T ML . Obviously, T Ml _ is equal 
to sea surface temperature because of the vertical homoge- 
neity assumption in the mixed layer. 

The surface boundary condition on heat flux is given 
by: 

(w'T 7 ) ^ - — — (2.6) 

f°« C P 



where the zero subscript indicates the value of w' T' at z=0. 
The net surface heat flux, Q/^Cj,, is prescribed to be 
harmonic in time : 



Q 





5^ " 



( 



■ 2 . rrv 

i 'Y c 



(2.7) 



The amplitude of this sine function is set equal to a 
representative mid-latitude value of 4’-10 -3< ’C cm/sec. In 
this study, heat flux is positive downward. Time t=0 is the 
start of the "heating season", when Q first becomes 
positive . 

The boundary condition on heat flux at the base of 
the mixed layer is given by the entrainment heat flux. 
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( 2 . 8 ) 



( w'T' ) K - ~ AT Vv' e , 

where the subscript, h, indicates the value of w 7 T 7 at z=-h, 
the base of the mixed layer. The variable. We is the 
entrainment rate and & T is the temperature jump across the 
entrainment zone at the base of f hi e mixed layer. The turbu- 
lent heat flux is linear with depth in the mixed layer (due 
to the assumption of a well-mixed layer), mcnotonically 
decreases with depth in the entrainment zone and is zero 
below. The profile of turbulent heat flux is shown in 
Figure 2.1. The net surface heat flux varies from positive 




Figure 2.1 Turbulent Heat Flux Profile. 

to negative values according to equation (2.7). The 
entrainment heat flux term, -iiT We, is always negative or 
zero. This is because the entrainment rats is positive when 
the mixed layer is deepening but "turns off" and is zero for 
non- entrain ing shallowing cases. The thickness of the 
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entrainment zone, 6 , is assumed to be much less than the 
depth of the mixed layer, h. 

Equation (2.5) now becomes 

^ ^rn. __ G AT We (2.9) 

a T /% H K 



Seta that the advecticn term is absent. Thus, 
motion does not directly affect the mixed layer 
Instead, the interior motion affects the 
temperature indirectly through the variables h. 



the interior 
temperature, 
mixed layer 
AT and We. 



2: C 1 



-IS - 



-aw 







Figure 2.2 Initial Temperature Profile. 
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The initial temperature profile is given in Figure 
2.2. The reference depth, D, is 200 m and will be used in 
some of the equations ir. later sections. 

2. Scale A na lysis and Simplification of t he Heat 

Equ ati on 

Consider now the relative magnitudes of the terms of 
equation (2.4) . In the mixed layer, the temperature may 
cnar.ge by as much as 10 °C in a half year. The time rate of 
change of mixed layer temperature would then be or. the order 
of 1 0 " 6 ° C/sec or less. Below the mixed layer the time rate 
of change of temperature of fluid particles will generally 
be smaller. 

The first term on the right of equation (2.4) is the 
turbulent heat flux term. In the mixed layer, the only term 
on the right side of equation (2.5) is the turbulent heat 
flux. Thus this term will scale" with the seasonal time rate 
of change of mixed layer temperature, 10 _6 °C/sec. The 
prescribed surface hsat flux boundary condition, Q/ p c c p , has 
extreme values of t4‘10" 3 °C cm/sec. The value of h varies 
from around 30 m to around 100 m. Thus, the first term on 
the right side of equation (2.9) has a range of values 
between ±10 -6c C/sec. The second term on the right of equa- 
tion (2.9) has roughly the same order of magnitude, but it 
can only be negative. Below the mixed layer and entrainment 
zone, that is, below the flux region, the turbulent heat 
flux ten vanishes. 

Next consider the advection term in equation (2.4). 
Recall that this term is zero in the mixed layer. The 
gradient of the temperature profile just below the mixed 
layer is initially about 10 _3 °C/cm. During the heating 
season the temperature at the top of the thermocline will 
increase, but the vertical extent of the thermocline will 
also increase as the mixed layer shallows. This term should 
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thus be expected to have the same order of magnitude. The 
amplitude of W will be prescribed according to equation 
(2.19) in section D. The maximum value at 200 m depends on 
the wave or oscillation frequency, but it will be on the 
order cf 10 -3 cm/sec. This amplitude will be smaller at 
lesser depths. The result is that the advection term of 
equation (2.4) should be of the order of 10 -8 °C/sec below 
the mixed layer. 

Finally, consider the diffusion term of equation 
(2.4). This term vanishes in the mixed layer. Below the 
mixed layer a reasonable scale for the second derivative of 
temperature is on the order of 4 c C/(5000 cm) 2 . A typical 
range of values for the eddy diffusivity, K, is 0.1 to 
0.5 cm 2 /sec. Thus, the diffusion term has a magnitude of 
about 10~ 8 to 10 -7 °C/sec. The diffusion term 'is much 
smaller than the other terms of equation (2.4) and will, 
from now on, be neglected. The heat equation (2.4) thus can 
be simplified to 




C. MIXED LAYER DEPTH 

In the model of Garwood (1977), the mixed layer either 
shallows to a retreat depth, h r , or deepens at the entrain- 
ment rate. We. The retreat depth is determined by assuming 
steady state conditions for the turbulent kinetic energy 
budget, letting entrainment heat flux vanish and letting 
<w" /2 > approach zero. The resulting equation for retreat 
depth is 
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= 



/ _ \ ^ / i. 

\ £ / 



( 2 . 11 ) 



3 /-*• ) C^)o 

(Garwood, 1977). In this equation, m , and m } are dimension- 
less model constants, <* is the thermal expansion coefficient 
and <E> is the mixed layer mean turbulent kinetic energy. 
Also, u„, the friction velocity, is the square root of the 
surface wind stress divided by density. Throughout this 
study, it is set equal to one centimeter per second. Since 
mixed layer mean turbulent kinetic energy depends on surface 
boundary conditions, the retreat depth is independent of the 
prescribed interior motion. 

In the case of active entrainment, the mixed layer deep- 
ening rate in the model depends on the entrainment rate and 
the mean vertical velocity at the bottom of the mixed 
layer due to the prescribed interior motion: 



ih_ 

J t 



: W, - W, 



( 2 . 12 ) 



The entrainment rate. We, depends on the entrainment zone 
temperature jump, AT, and the entrainment heat flux (see 
equation (2.8)). The entrainment heat flux is modelled by 



(?rV ^ 

h 



( 



1/1 / S' \ 
/ 



'M Nt 



(2. 13) 



*3 



(Garwcod, 1977), where 



is a model constant. 



The 



entrainment rate can then be expressed as 



+ m 4 ( w ' i 



‘/a. / - X 

<E > 



w. 



*3 



U AT 



(2. 14) 
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Therefore, equations (2.12) and (2.14) give a solution for 
dh/dt if W^, <vT r 2> and <E> are known. 

For moderately strong wind stress, mixed layer mean 
turbulent kinetic energy, <!>, depends primarily on surface 
wind stress, which is held constant here. The vertical 
component of turbulent kinetic energy, (w T ”2) 1/2 , depends 
primarily cn surface heat flux and surface wind stress. 
When the surface heat flux becomes large enough, this term 
becomes zero, stopping the entrainment process and permit- 
ting shallowing. These two terms are unaffected by the 
interior motion. The entrainment rate is affected by the 
mean vertical motion indirectly as a result of changes in 
mixed layer depth, h, and entrainment zone temperature jump, 
AT. 



D. I NT EH 10 R VERTICAL MOTION 



Consider for the moment a two layer, deep ocean with the 
•upper layer about 200 m deep and the lower layer about 
5000 m deep. Assuming the fluid is incompressible, the 
continuity equation is differentiated with respect to z: 



6? 



w 






u 



- £L_ v , 



(2. 15) 



where U, V and w are the mean velocity components in either 
layer. Furthermore, assume that horizontal momentum is 
vertically well mixed in the upper layer. Thus, 

c) u J77 

— =■ O ; - 0 (2. 16) 
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w it h a 



Than N must ba linear with depth in tne upper layer, 
rigid-lid approximation, W = 0 at the surface. 

Our interest is confined to the upper 200 ra and a linear 
vertical profile for W is a good approximation for long 
period interior motion. The functional dependence of w on 2 
is 

W (O - - a t (2.17) 



where D = 200 m is the 
interior vertical motion 
function of time, namely 



reference depth and is the 

at 200 m given by a sinusiodal 



w d - 



- W 0 si 



/ iTTt 
l“T 




(2. 18) 



The amplitude, W 0 , is defined by 



w 



O 



•ITT D \ O Kva 

T 1 O rr> 



(2. 19) 



For the one-year period interior motion, 

W 0 = 8*10~* cm/sec - 69 cm/day. The value of W c will be 2, 
3 and 4 times this for the one-half, one-third and one- 
guarter year periods, respectively. These values of W 0 are 
selected so that the displacement amplitude at 75 m is 
approximately 15 m. The actual values of the displacement 
amplitude depend on phase and will be discussed below. 

Combining equations (2.17) and (2.18) yields the 
equation for the prescribed interior motion. 



w ( * , ± ) 



a 






w„ 







( 2 . 20 ) 
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The formula for the particle position due to the interior 
motion can be obtained by noting that W = dz/dt, and is 
given by 



i 






-TW, 
a-rr D 




( 2 . 21 ) 



For any fixed phase the maximum and minimum values of z are 
given by 
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and 
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(2.23) 



The peak-to-peak displacement amplitude is the difference. 
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(2. 24) 



The maximum peak-to-peak displacement amplitude will occur 
when cos p =+1 and the minimum when cos <p =- 1 . For the 
vertical motion amplitude given by equation (2.19) , a 
particle with an initial position at z=-75 m will experience 
a displacement amplitude (half of the peak-to-peak value) 
between 12.4 and 18.4 a. 

The model time step is one hour. This is much smaller 
than the cne-year or even one-quarter year periods consid- 
ered. Hence, since W ^ does not vary significantly over a 
one hour time step, the integral of (2.17) can be approxi- 
mated by 
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where /}t=t t -t, is the time increment. (one hour) and is 



obtained from equation (2.18). In the model, vertical 
advection is handled by compressing or expanding the 
vertical grid according to equation (2.25). No significant 
difference in the results was found between model runs using 
equation (2.25) and equation (2.21). 

The vertical increment is initially one meter, but this 
varies as a result of the compression or expansion of the 
vertical grid. This variation in the vertical increment is 
never more than about 25^. A consequence of this scheme is 
that the lowest grid point is not always at a depth of 199 m 
as in the initial profile. Figure 2.2. Care must be taken, 
therefore, when performing vertical integration from the 
surface to the reference depth, D = 200 ra. The integrals 
that will be calculated in the next section involve the 
temperature difference, T -T 0 . Problems with varying the 
grid size can easily be avoided by setting the reference 
temperature, T 0 , equal to Tx>, the initial temperature at 
200 m. This way, T-T D will always be zero at the base of 
the vertical grid. The grid is initially made deep enough 

so that the seasonal thermocline never reaches the bottom of 
the grid; hence, f D always equals its initial value, 6.08°C. 

Subtracting the reference temperature, f D , from tempera- 
ture, T, in the integrals in the next section amounts to 
subtracting the terms: 
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and 
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Since f p always remains constant, both expressions, (2.26) 
and (2.27), will always remain constant. Subtracting T^ 
from T will amount to subtracting constants from the 
integrals if the next section. 

E. POTENTIAL ENERGY AND HEAT 

The potential energy per unit area for a column of water 
from the surface to depth, D, is given by 
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(2 . 28) 



Employing equation (2.2) with as the reference 

temperature, allows equation (2.28) to be written as 
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The constant in equation (2.29) is 
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(2. 30) 



The integral in equation (2.29) will be used for the compar- 
ison of the potential energy with and without the interior 
motion. It will be represented by the variable P with units 
of 0 C m 2 where 



P = - ( ( T - Tj, ) 2 cU . 



(2. 31) 
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Similarly , 
with units 



heat content will be indicated- by the variable H 
°C o' where 
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(2.32) 



The analysis of the physical effects of the interior 
motion can be facilitated by developing the expressions for 
the time rate of change of H and of P. First, 
differentiating equation (2.32) with respect to time gives 
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Equation (2.33) can be written as 
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Substituting equation (2.10) for the time rate of change of 
temperature into equation (2.34) yields 
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(2. 35) 



The first term on the right , when integrated, is just the 
net surface heat flux. Substituting equation (2.17) into 
the second term on the right of (2.35) and integrating by 
parts gives 
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(2. 36) 
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Using equation (2.26) this can be rewritten as 
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Equation (2.35) can new be written as 
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The formulation for the rime rare of change of the 
potential energy is similar and is given by 
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Substituting equation (2.10) into (2.39) yields 
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(2. 40) 
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(2. 41) 
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The first term on the right of equation (2.41) vanishes 



because w'T' = 0 at z = -D. Since w'T' is zero below the 
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entrainment zone, 



equation (2.41) reduces to 
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This integral may be broken up over the mixed layer and the 
entrainment zone, producing 
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where < w ' T 7 > is the vertically averaged mixed layer mean of 
the turbulent heat flux. Since w ; T' is linear with respect 
to z in the mixed layer, this mixed layer mean can be 
written 
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Consider* t h e second term on the right side of equation 
(2.40). Substituting equation (2.17) for W and integrating 
by parts yields 
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using equation (2.27) this can be rewritten as 
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Cn combining equations (2. 4 0), (2.43), (2.44) and (2.46), 

the time rate of change of the potential energy variable now 
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Because the turbulent heat flux is monotonic across the 
entrainment zone, as illustrated in Figure 2. 1: 



This indicates that the integral term in equation (2.47) 
should generally be smaller in magnitude than the first term 



entrainment heat fluxes nearly cancel. 

Actual values of the terms of equation (2.47) were 
calculated for one -year period interior motion at two 
different phases and the second term on the right of the 
equation was smaller than the first, usually by one order of 
magnitude. Furthermore, the third term on the right was 
generally an order of magnitude larger than either of the 
other two. Equation (2.47) will thus be approximated by 
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on the right side, unless of course, the surface and 
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III. EXPERT M ENTS AND RESULTS 



A. SUMMARY OF PERTINENT EQUATIONS 

The equations for time rate of change of mixed layer 
depth, mixed layer (and sea surface) remperarure , heat and 
potential energy are 
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The retreat depth is given by 
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where (w'T') has been replaced by equation (2.6) 
The entrainment rate is given by 
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The prescribed net surface heat flux is given by 
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and is plotted as a function of time in Figure 3.1. 




The prescribed vertical motion at the reference depth, 
D=2 0 0 m, is 
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where W 0 is defined by equation (2. 19) . 
z=-h into equation (2.17), it can be seen 
directly as h and W D . 

Some additional notation will be useful 
sion of the results. To observe the effect 



3y substituting 
that W K varies 

for the discus- 
of the interior 
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notion, the values of h, T fv , L , H and ? are compared to 
without interior motion. Contours are plotted of 
difference of mixed layer depth with and without 
interior motion, likewise for mixed layer temperature, 
and potential energy. Therefore, the following 
variables are defined: 

A K = U - k I (3 

1 W £ 0 ’ 



AT 



ML 






ML 
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A P - ? 



z 0 
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B. RESULTS WITH HO INTERIOR MOTION 

The mixed layer depth has been plotted in Figure 
Note that the vertical axis scale is inverted for cla 
The year was started at the beginning of the heating s 
and the mixed layer retreats for the first quarter 
Then the mixed layer deepens slowly for the next half 
and more rapidly in the last quarter. In the last few 
the mixed layer shallows to the t=0, 1 year,...,n- 

cyclical steady state value of about 75 m. 
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The mixed layer temperature, f , is plotted in Figure 
A. 5. As expected, the mixed layer temperature increases 
during the heating season and decreases during the coding 
season, for the most part. However, note that the peak 
mixed layer temperature occurs around day 170 and T H1 _ starts 
to decrease ten days before the end of the heating season. 
This is because the entrainment heat flux term has begun to 
dominate eguation (3.2). Some entrainment had already 
started fifty days earlier. By day 170, towards the end of 
the heating season, AT has reached its largest value, 
7.8 °C . Cn day 170 the cooling of the mixed layer due to 
entrainment has reached 0. 021 °C per day while the surface 
heating has dropped to 0.030 °C per day. 

The heat content, plotted in Figure A. 6, depends only on 
surface heat flux for the case of no interior motion (see 
eguation (3.3)). Any heat lost from the mixed layer due to 
entrainment will be gained by the region below. Since heat 
is calculated from the surface to 200 m, entrainment will 
have no effect on its value. 

Potential energy per unit area is plotted in Figure A. 7. 
This, too, generally increases with heating and decreases 
with cooling. However, note that the peak value occurs at 

about day 305, 25 days into the cooling season. This is 

because potential energy is also increased by mixing as 
indicated by eguation (3.4) . The entrainment term dominates 
eguation (3.4) from day 170 until day 305, when the net 
surface heat flux is near zero. A close inspection of 
Figure A. 7 reveals that the potential energy at the end of 
the year is greater than the initial value. Specifically, 
the initial value of F is 19608 Z 'C m 2 ; the value at the end 
of the first model year is 20204 °C m 2 and at the end of the 
second year is 20579 ’C m 2 . This increase in the potential 
energy is due to the adjustment cf the temperature profile 
as the ocean approaches cyclical steady stats. 
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Figure A. 1 is a plot cf temperature versus depth and 
time. The top of the graph where the temperature contours 
are vertical is the mixed layer. The mixed layer structure 
is quits evident after the first quarter year due to the 
strong temperature gradient at the base of the mixed layer. 
Note that the temperature profile at any depth below the 
mixed layer is left unchanged until mixing reaches that 
depth . 

C. ONE-YEAR PERIOD INTERIOR MOTION 

The zero contours of the prescribed interior motion have 
teen drawn in Figure A. 8 with upwelling and downwelling 
indicated. The dotted contours indicate when the magnitude 
cf the interior motion drops to half of the maximum value. 

During the first quarter year the mixed layer is shal- 
lowing. Since the retreat depth depends primarily on 
surface boundary conditions, mixed layer depth will be unaf- 
fected by interior motion at this time. This is borne out 
by the contours of /ih (defined in equation (3.9)), Figure 
A. 10. The first quarter year shows no difference in mixed 
layer depth. With little or no effect on mixed layer depth 
during the first quarter year, there will be little cr no 
difference in the value of the first term of equation (3.2) 
with and without the interior motion. The second term of 
equation (3.2) will be zero at this time since entrainment 
stops during shallowing. Hence, mixed layer temperature 
will also remain unaffected by the interior motion during 
the first quarter year. This can be seen in Figure A. 12. 
Heat and potential energy, on the other hand, are influenced 

by the advection of temperature in the thermocline as indi- 
\ 

cated by the presence of the last term in equation (3.3) and 
equation (3.4). Thus, the interior motion will influence 
heat content and potential energy even during the first 
quarter year (Figures A. 16 and A. 18). 
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1 • Eff ect on M i xed Lay.er De p~.h 

Equation (3.1) indicates •chat downwelling will 
enhance deepening, and upwelling will reduce deepening. In 
Figure A. 9, it can be seen "hat deepening still dees occur 
after the first quarter year as it did without any interior 
motion. However, the rate of deepening depends on phase. 

At the 90° phase, there is upwelling during "he 
middle half of the model year (see Figure A. 2) reducing the 
rate of deepening. At day 274 down welling begins and the 
mixed layer begins to deepen more rapidly. The very rapid 
deepening in the last quarter of the year for this phase 
results from both entrainment and downwelling working 
together. Calculated values of the terms of equation (3.1) 
indicate that entrainment is the dominant term most of the 
time during deepening. In this example the entrainment rate 
has bean enhanced by conditions set-up by the interior 
motion. Equation (3.6) indicates that entrainment rate. We, 
is inversely proportional to mixed layer depth, h, and 
entrainment zone temperature jump, at. Compare Figure A. 2 
to Figure A.1 (the case with no interior motion) and A. 3 
(the 270° phase case). It can be seen that at the start of 
the last quarter of the yaar the 90° phase case has a much 
lower mixed layer temperature than the other two cases (the 

reason for this will be discussed in the next section) . 

This lower mixed layer temperature means that the entrain- 
ment zone temperature jump will be smaller, allowing for a 
larger entrainment rate. At the start of the last quarter 
of the year, the differences in h are not as significant as 
those in AT. 

For example, the values of at, h. We and W k have 

teen calculated for these three cases for day 280. The 

mixed layer depth, h, is about 50 m for the case without 
interior motion, 65 m for the 90 c phase and 68 m for the 
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270° phase. These values dc not greatly differ. However, 
the temperature jump for the 90' phase case is 1.3°C which 
is less than half the 2.9°C value for the case with no 
interior motion which in turn is less than the 4.7 c C value 
for xhe 270° case. The entrainment rats for the 90 1 -' phase 
case is 90 cm/day compared to 43 cm/day for xhe no interior 
motion case and 22 cm/day for the 270° case. For both 
phases considered, the interior motion is at least an order 
of magnixude smaller than the entrainment rate. 
Nevertheless, at day 280, the deepening in the 90" phase 
case is helped a little by downweliing and in xhe 270° case, 
hindered a little by upweiling. 

In the last few days of the year, xhe increase in 
mixed layer depth begins xo take over, particularly for the 
90° case. This inhibits further entrainment and the mixed 
layer then rapidly shallows to the cyclical sxeady-sxate 
value of 75 m (Figure A. 9) . 

In the 270 ° phase case (Figure A. 3) there is down- 
welling during the middle half of xhe year, which increases 

the deepening of the mixed layer. By the sxart of the last 
quarter of xhe year, the mixed layer is warmer than in the 
case with no interior mcxion or in the 90 c phase case. This 
leads to a larger value of A T and a reduced entrainment 
rats . 

Comparing Figures A. 8 and A. 9 it can be seen x hat 
mid-year upweiling is associated wixh a reduced raxe of 
deepening and mid-year downweliing is associaxed with an 
increased deepening rate. Looking at Figure A. 10, it can be 
seen that the actual decreased or increased values of mixed 
layer depth lag the start of upweiling or downweliing 
(respectively) by 90°. The time rate of change of mixed 
layer depth tracks along with the interior motion in 
accordance with equation (3.1). However, acxual mixed layer 
depth is the integral of this and integrals of periodic 
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functions lag them by 90° . Also note -hat the strong 
enhancement of and strong hindrance to deepening at the end 
of the year (Figure A. 10) lags the strong decrease and 
strong increase in mixed layer temperature (Figure A. 12) by 
about 90 c . 

2. Effect on nixed Layer T emper a ture 

Consider now the mixed layer heat equation (3.2). 
The net surface heat flux is prescrioed and unaffected by 
the interior motion. The entrainment heat flux, given by 
equation (3.6), is inversely proportional to mixed layer 
depth. The other terms in equation (3.6) depend primarily 
on surface boundary conditions. So, the first term in equa- 
tion (3.2) is inversely proportional to mixed layer depth, 
h, and the second, inversely proportional to the square of 
h. This implies that an increase in h will reduce the 
magnitude of mixed layer temperature change; a decrease in h 
will increase the change in T^,_. This is because, with a 
deeper mixed layer, a given amount of heat added to or 
removed from the mixed layer must be spread over a larger 
mass cf water. 

This effect is illustrated in the examples shown in 
Figures A. 2 and A. 3. In Figure A. 2, the mixed layer is 
shallower than in Figure A. 3 from about day 180 to about day 
300. This resulted from the mid-year upwelling for the 90 s 
phase case and the aid-yea r downwalling for the 270° case. 
The shallower mixed layer in Figure A. 2 allows for a much 
more rapid temperature drop than occurs with the deeper 
mixed layer of Figure A. 3. This can be seen in the enhanced 
gradients in Figure A. 2. After day 300, the rapid deepening 
of the mixed layer for the 90° phase case and slew deepening 
for the 270° case reverse the situation. 

To see the effect of the interior motion for all of 
the phases, the T.vu. contours should not be compared directly 
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to the interior motion contours but rather to the contours 
of the mixed layer depth difference, Ih (defined by equa- 
tion (3.9)). Between the 0 C and 133° phases, after day 180 
there is a region in Figure A. 10 where the mixed layer is 
shallower with the interior motion than without." This 
causes the mixed layer temperature decrease (the second half 
of the year is the cooling season) to be greater, resulting 
in the negative values of /_iT ML in the last quarter year at 
these phases (Figure A. 12). Between the 180° and 360 ? 
phases the opposite occurs. 

At the end of the year the mixed layer depth returns 
to 75 m. Mixed layer temperature does not return to its 
initial value. The second year T^ L data is plotted in 
Figures A. 13 and A. 14. The reduced temperatures for phases 
0° to 190" and increased temperatures for phases 180° to 
360° are continued into the second year. The mixed layer 
depth follows the same pattern in the second year as it did 
in the first. This causes the f hL pattern begun in the 
first year to be enhanced in the second year. Continued 
pumping of the ocean by this period interior motion enhances 
this structure. If permitted, it would eventually cause the 
upper 200 m to become isothermal (in the region of the 90° 
phase at the end of the year) , followed by deepening beyond 
the limits of what the model can handle. It would not be 
realistic, however, to have this single period interior 
motion continue for the several years necessary for this to 
happen. In reality, the interior motion would probably be 
damped by interactions with the boundary layer mixing 
processes. There should also be more than one period 
present. Furthermore, surface heat flux should be influ- 

enced by sea surface temperature. These three possibilities 
open up avenues for further research. 
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. E if set on Hear Content 

Equation (3.3) mdica-.es that. change in the heat 
content depends on net surface heat flux and on pumping of 
the ocean by the interior motion. The net surface heat flux 
is prescribed and is independent of the interior motion. 
Subtracting equation (3.3) for the case of no interior 
motion from the same equation with interior motion gives 

<•-' ^ ^ _ _ w D (3. 13) 



where AH is defined by equation (3.11). Equation (3.13) 

indicates that upwelling will cause AH to decrease and 
downwelling will cause AH to increase. This can be 
observed by comparing Figures A. 8 and A. 16. The input to 
heat content due to net surface heat flux alone causes H to 
increase during the first half of the year and decrease 
during the second half as shown in Figure A. 6. The combina- 
tion of surface heating with the effect of the interior 
motion produces the saddle-shaped curve shewn in Figure 
A. 15. 

The effect of different phases of the interior 
motion can be demonstrated by the following two simple exam- 
ples (Figures 3.2 and 3.3) . No mixing is permitted and a 
two-layer ocean is assumed for simplicity. 

In the first example shown in Figure 3.2 there is 
initial downwelling followed by heating, then upwelling 
followed by cooling and final downwelling to the initial 
mixed layer depth. When downwelling precedes heating, the 
heat is distributed ever a thicker mixed layer. A fixed 
amount of heat will produce a smaller increase in T than 
would have resulted without the downwelling. If this is 
followed by upwelling before cooling, the thinner mixed 
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Figure 3.2 Heat Content Sxanple 1. 
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Figure 3.3 Heat Content Example 2. 
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layer will reach a lower temperature for a given amount of 
heat removed. At the end of the complete cycle the mixed 
layer temperature will be lower. Recall that the heat 
content is defined by 

r° 

hi - ) ( T - T & ) a i . 

~ D 

Therefore the lower mixed layer temperature produces a lower 
heat content. 

Upwelling preceding the heating season distributes 
the heat ever a thinner mixed layer. This allows a greater 
increase in P than would otherwise occur. If this is 

followed by downwelling preceding the cooling season, a 
fixed amount of heat removed will produce a smaller drop in 
T. The net effect will be a higher temperature in the mixed 
layer and a higher heat content. 

Because of the asymmetry, as illustrated in these 
examples, heat is net conserved for a particular phase of 
the interior motion. This can be seen in Figures A. 15 and 
A. 16. At some phases the year-end heat content is higher 
than it was at the beginning of the year, at other phases, 
lower. However, if the heat content is averaged over all of 
the phases, the year-end value equals the initial value. 
The average heat content for sixteen equally spaced phases 
from 0 C to 360° was tabulated. The phase-averaged value of 
H was 0.07 C C m less at the end of the year than it was 
initially. This small error can be attributed to the coarse 
resolution of averaging only sixteen phases over the entire 
360 7 . 

4 • E f f jet 2.2 Pote nt la l Energy 

Time rate of change of potential energy is given by 
equation (3.4). From this equation it can be seen that 
surface heating tends to increase potential energy, surface 
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cooling tends no decrease it. The net surface heat flux is 
independent of the interior motion so the effect on iP due 
to this term depends on <clh. If there is an increase in 
mixed layer depth during heating, the buoyancy added will be 
mixed deeper, increasing potential energy. If there is an 
increase in mixed layer depth during cooling the effect will 
be opposite. 

The entrainment heat flux term has a positive sign 
in equation (3.4). Increased mixing moves buoyant water 
down and less buoyant water upward, increasing potential 
energy. 

Comparison of values of the terms of equation (3.4) 
indicates that the third term on the right is generally 
dominant. The exception is when the interior motion goes to 
zero, twice a year. The first term of equation (3.4) will 
not differ greatly with and without interior motion unless 
there is- a large difference in mixed layer depth. Recall 
that entrainment heat flux is nearly inversely proportional 
to mixed layer depth (equation (3.6)). Thus, the second 
term will net vary significantly. So, the potential energy 
closely follows the interior motion. This can be seen by- 
comparing Figures A. 8 and A. 18. During downwelling, 
potential energy increases, during upwelling it decreases. 

The physical reasoning behind the effect of the 
third term of equation (3.4), the pumping term, is that 

downwelling pushes buoyant surface water downward, whereas, 
upwelling is in the same direction as the buoyant forces. 

5. Summary of Respo nse to Period In t er io r 

A ction 

During the first quarter year there is shallowing 
and the interior motion has little effect on either mixed 
layer depth or temperature. Mid-year upwelling reduces 
deepening and mid-year downwelling enhances deepening. The 
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shallower mixed layer resulting from upwelling will cool 
faster; the deeper mixed layer resulting from downwelling 
will cool slower. A cooler mixed layer has a reduced jT 
allowing for a larger entrainment rate and more rapid 
deepening. A warmer mixed layer deepens more slowly. 

Potential energy and heat are both directly affected 
by the interior motion. Beth are increased by downwelling 
and decreased by upwelling. 

Before examining the mixing response tc higher 
frequencies, an observation can be made regarding horizontal 
structure produced by the interior motion. As noted in 

chapter one, phase may be related to horizontal displacement 
by substituting 
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(3. 14) 



into equation (2.20) to obtain 
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(3. 15) 



The wavelength, L, might be on the order of 300 km, large 
enough for the one-dimensional assumption to be valid. 

Figures A. 1 2 and A. 14 show the development of hori- 
zontal bands of positive and negative mixed layer tempera- 
ture anomalies- which would result from a one-year period 
interior wave. Because temperature is assumed to be verti- 
cally homogeneous in the mixed layer, T Ml also represents 
sea surface temperature. Figures A. 11 and A. 13 show a time 
sequence of the horizontal sea surface temperature field 
resulting from the interior motion, the heat cycle and 
mixing processes. The entire surface warms during the 
heating season and cools during the cooling season. During 
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the first half of the first year, there is no significant 
difference in sea surface temperature along the horizontal 
axis. After the first half year, the region where there is 
upwelling at the start of the cooling season becomes and 
remains cooler than the region where there is dcwnwelling at 
the start of the cooling season. 



D. EFFECT OF SHOBTEB PEBIOD INTEBIOB MOTIOH 

Comparison of Figures A. 19 and A. 21 shows that as before 
downwelling increases deepening and upwelling reduces deep- 
ening. This can also be observed by comparing Figures A. 30 
and A. 32 and by comparing Figures A. 39 and A. 41. The mere 
rapid variation in mixed layer depth dees not yield the 
strong effect on temperature found for the one-year period 
interior motion. This can be seen in the weaker gradients 
of £T h( _ cf Figures A. 23, A. 34 and A. 43 compared to that of 
A. 12. Hence, the higher frequency interior motion does not 
have as strong an effect on the mixed layer as the lower 
frequency metien does. 

The second year temperature data was plotted for the 
half-year period interior motion in Figures A. 24 end A. 25. 
Notice that the same structure seen in the first year reap- 
pears in the second, but the cold spots have not grewn 
significantly colder nor the warm spots warmer as in the 
case of the one-year period motion. 

As before, the heat and potential energy differences, 
AH and AP, track along with the interior motion as seen in 
Figures A. 27, A. 29, A. 36, A. 33, A. 45 and A. 47. The 
resulting contours of H and P appear more complex. This is 
because the values cf H and P without the interior motion 
follow a one year cycle (Figures A. 6 and A. 7) and the higher 
frequency cycles due to the pumping by the interior metion 
are superimposed on the one year cycle. 
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Fcr the cne- yea r period interior motion, upwelling lasts 
long enough'to enhance cooling significantly (at phases 
around 90°) and downwelling lasts long enough to inhibit 
cooling significantly (at phases around 270°). The shorter 
periods dc not maintain upwelling and downwelling long 
enough no produce an effect this strong. 

E. ASYMMETRY IN ENHANCED AND REDUCED MIXING 

In Figures A. 48 and A .49 the values of peak enhanced 
cooling and peak reduced mixed layer cooling are plotted. 
Both first and second year values are considered (in Figures 
A. 48 and A. 49, respectively). Negative values of a T iU 
indicate enhanced coding due to the interior motion. The 
maximum absolute value of all of the negative a'T mi . data 
points was determined to indicate the peak enhanced cooling. 
Positive values of at m1 _ indicate reduced cooling. The 
maximum value of all of the positive aT wl data points indi- 
cates the peak reduced cooling. The positive values dc not 
generally indicate enhanced heating, because the wave has 
very little effect on the mixed layer temperature during the 
heating season. 

The first obvious observation that can be made from 
Figures A. 48 and A. 49 is that the largest peak anomalies 
occur with the one-year period interior motion. For the 
cne-year period interior motion the peak enhanced and peak 
reduced mixed layer cooling is 2.04 and 2. 16°C, respec- 
tively. For the one-half year period motion these values 
are 1.35 and 1.15°C, respectively. Also, for each of the 
periods studied, the difference between the peak enhanced 
cooling and the peak reduced cooling is greater in the 
second year than in the first. All of the periods studied 
except the one-year period have greater peak enhanced 
cooling than peak reduced cooling. 
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The phase-averaged values of rhe variables Ah, AI X , 
AH, ana A? were computed for sixteen equally spaced phases 
from 0° to 360° . The phase-averaged values at the end of 
the first and the end of the second model years are listed 
in Table I. These values indicate the net effect of the 
interior mc-ion after each year of interacting with the 
mixing processes. Consistent with the resulrs shown in 
Figures A. 48 and A. 49, the results in Table I indicate that 
the one-year period interior motion has produced a net 
reduction of mixed layer cooling. The higher frequency 
interior motion has produced a net enhancement of mixed 
layer cooling. For all of the periods studied there has 
been a net reduction of the mixed layer depth. Heat is 
conserved, hence the phase -averaged value of AH is zero at 
the end of each year. The effect of the one-year period 
interior motion was to decrease the potential energy. The 
higher frequency interior motion produced a net increase in 
potential energy. 

A net - increase (decrease) in potential energy indicates 
that the interior motion has enhanced (reduced) the overall 
mixing. The amount of deepening required to produce an 
equivalent increase in potential energy is instructive. 
This can easily be calculated for the simplified example of 
a two-layer ocean, illustrated in Figure 3.4 with initial 
mixed-layer temperature, T 0 and depth, h a . This mixed layer 
is deepened to h, with new mixed layer temperature, T, . 
Heat is conserved, so, using the definition for the heat 
variable, H, given in equation (2.32) : 

( 1 i “ 1 n ) K = ( ' o - 1 t> ^ (3. 16) 

Using equation (2.31), the change in the potential energy is 
given by 
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(3.17) 



P, - ? 0 = ( T, - T* ) C _ C T„ - T 0 ) UT . 

a • i 

Sabstituting equation (3.16) into aquation (3.17) yields 

?, - K = ^(T c - T„) Uo( w, - k.) (3.18) 

The term T„-T D is just the initial entrainment zone tempera- 
ture jump. The increase in potential energy due to deep- 
ening with boat conserved is directly proportional to the 
increase in mixed layer depth for the two-layer example. 



— > T 





initial profile after mixing 

Figure 3.4 Two-Layer Ocean Example. 

Consider now the pha se- averaged potential energy anoma- 
lies at the end of the first year, shown in Table I. 
Suppose the initial mixed layer depth is 75 m and the 
initial temperature jump is 4 °C. The net reduction of 
potential energy for the one-year period interior motion 
would then be equivalent tea reduced mixed layer depth of 
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0.73 m. The net increase of potential energy for the or.o- 
half, one-third, and one-quarter year period interior motion 
is equivalent to an increase of mixed layer depth by 1.6, 
1.8 and 1.8 m, respectively. 

The amount of dcwnweiling and upwelling are the same 
over the whole year and for all phases. Yet, the net reduc- 
tion of mixing and mixed layer cooling does not equal the 
net enhancement of the same. The net effect of the one-year 
period interior motion is to hinder mixing. The net effect 
of the shorter period interior motion is to enhance mixing. 



TABLE I 

Phase-Averaged Hean Anomalies 





Phase-averaged mean of: 




Da y 


Ah 




AH 


AP 


Period 




(m) 


( C C) 


(*C m) 


rc m ) 


(years) 


0 


0. 


0. 


0. 


0. 




36 5 


-0. 43 


+ 0. C6 


0. 


- 109. 


1 


73 0 


-0.71 


+ 0. 22 


0. 


-645. 




0 


0. 


0. 


0. 


0. 




365 


-0. 16 


-0.05 


0. 


239. 


1/2 


730 


-0. 16 


-0. 06 


0. 


33 1 . 




0 


0. 


0. 


0. 


0. 




365 


-0. 23 


-0. C5 


0. 


26 7. 


1/3 


730 


-0. 23 


-0. 07 


0. 


435. 




0 


0. 


0. 


0. 


0. 




365 


-0. 31 


-0. 05 


0. 


27 5. 


1/4 


730 


-0. 31 


1 

O 

« 

o 


0. 

J 


451. 
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IV. CONCLUSIONS AND R ECOM fl EN DATIO NS 

A. CONCLUSIONS 

Seasonal-scale interior motion has a significant effect 
on boundary layer mixing processes. This effect is most 
pronounced whan the period of the interior motion is one 
year. For this period, upwelling between the heating and 
cooling seasons can strongly enhance mixed layer cooling and 
subsequent deepening. Downwelling can strongly reduce the 
same. This effect may produce a significant, observable 
horizontal pattern in sea surface temperature. Although 
there is negligible effect on sea surface temperature and 
mixed layer depth during the beginning of the heating 
season, the interior motion dees affect the potential energy 
and heat content at this time. Throughout the year the 
potential energy and heat are directly related to the 
pumping by the interior motion. Upwelling decreases and 
downwelling increases both potential energy and heat 
content. 

Higher frequency interior motion does not have as strong 
an effect. The cycling between upwelling and downwelling 
during the cooling season prevents the strongly enhanced 
mixed layer cooling from occurring. 

The effect of the interior motion is asymmetric. 
Generally, the enhancement of mixing and mixed layer cooling 
exceeds the reduction of the same. The exception to this 
occurs for the one-year period interior motion. 

B. BECOBHENDATIONS 

One aspect of this research which departs from reality 
is in prescribing a single wave for the interior 
This could be changed quite easily by prescribing 
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motion . 



w » = £ w 




(4.1) 



The amplitudes, W^, and periods, T„ , could be based on the 
first few dominant frequency components observed in actual 
wave spectra. 

Another way that reality could be more closely 
approached is by allowing the net surface heat flux to vary 
depending cn the sea surface temperature. This should 
dampen the strongly enhanced or reduced cooling of the mixed 
layer temperature seen with the one-year period motion. 

An area open to further research is to use a finite 
depth ocean and an appropriate vertical wave structure. An 
appropriate wave structure for planetary waves is given by 
Willmott and Mysak (1980) . In their work the vertical 
motion is linear in the mixed layer: 



but, below the mixed layer the vertical motion is given by 



where N 2 is the 3runt-Vaisala frequency and c^ represents 
the speed of propagating long waves in the non-rotating case 
for mode number n. The solutions of equations (4.2) and 
(4.3) must be matched at the base of the mixed layer, z=-h, 
to ensure continuity cf the vertical velocity. The boundary 
conditions are rigid lid and no normal flow through the 
ocean bottom, so, w=0 at z=0 and at the ocean bottom. At 

each time step the mixed layer depth is recalculated and so 
the vertical wave structure problem would also have to be 
solved for the new w (z) . The combination of this sort of 



d W 



(4.2) 
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vertical wave structure with the Garwood model for boundary 
layer mixing processes could provide an interesting area for 
future study. 

In this study the interior motion was prescribed and 
acted on the mixed layer. The interior motion was not 
allowed to be changed by the mixed layer processes. If, for 
example, the wave motion enhances mixing, the increase in 
potential energy should be balanced by a decrease in wave 
energy. Although the frequencies considered here are much 
lower than those in the work by Bell (1978), the processes 
involved could be similar. 
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APPENDIX A 
FIGDBES 

Solid contours indicate positive values. 
Dashed contours indicate negative values. 
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Figure a. 5 Mixed Layer Temperature, No Interior Motion. 
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PE 10000 C-M-M 

2.0 2.2 2.4 2.6 2.8' 3.0 HEAT 100 C"M 




TIME DAYS 



Figure A. 6 Heat, No Interior Motion. 




Figure A.7 Potential Energy, No Interior Motion. 
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Figure A. 9 h (meters) , One-Tear Period Interior Motion. 




Figure A. 10 Ah (meters). One- Year period Interior Motion. 



60 




Figure A. 11 t iM1 _ (C) , One-Year Period Interior Hotion. 




Figure A. 12 AT Ml _ (C) , One-Year Period Interior Hotion. 
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Figure a. 15 H (lOO'C ■) , One-Year Period Interior Motion. 




Figure a. 16 AH (100 c C m) , One-Year Period Interior Motion. 
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Figure A. 17 P (10* °C m z ) , One-Year Period Interior Motion. 
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Figure A. 19 



1/2-Year Period Interior Hotion. 
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Figure A. 20 h (asters), 1/2-Year Period Interior Motion. 




Figure A. 21 



Ah (meters) , 1/2-Year 



6 6 



Period Interior Motion 






Figure A. 22 T Ml _ (C) , 1/2-Year Psriod Interior Motion. 




Figure A. 23 4 (C) , 1/2-Iear Period Interior Motion. 
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Figure A. 24 T Hl _ (C) , 1/2-Year Period, Second Year. 




Figure A. 25 AT, nt _ (C) , 1/2-Year Period, Second Year. 



68 






Figure A. 27 (100 V C a), 1/2-Year Period Interior Motion, 
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Figure A. 28 P (10* X a 2 ), 1/2-Tear Period Interior notion 



O 




Pigure A. 29 AP (10* °C m 2 ), 1/2-lear Period Interior Moti 
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Figure 4.30 1/3 -Ye ar Period Interior Motion. 
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Figure A. 31 h (asters) , 1/3-Year Period Interior Hotion. 




72 





PHASE DEG c PHASE DEG 

135 180 225 270 315 360 £ 9 < 5 9 ,° *? S l ?° 225 2? ,° 3 ! 5 3 ^° 





73 




Figure ft. 35 H (100 "C n) , 1/3-Year Period Interior Motion- 




Figure A. 36 



AH (100 C C m) 



1/3-Ysar Period interior 



Motion. 
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Figure A. 37 P (10* ‘C m 2 ) , 1/3-Tear Period Interior Motion. 




Figure A. 38 AP (10*°C a 2 ), 1/3-Tear Period Interior Motion. 
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Figure A. 39 



1/4-Year Period Interior Motion- 



76 





Figure A. 41 Ah («eters) , 1/4-Year Period Interior Motion, 
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Figure A.42 <C) , 1/4-Year Period Interior Motion 

O 




Figure A. 43 AT ML <C) . 1/4-Iear period interior Motion 
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PURSE DEG e PHASE DEG 




A. 44 H (1 00 ‘'C m) , 1/4-Tear Period Interior Motion. 




Figure A. 45 AH (100 c C m) , 1/4-Year Period interior Motion. 
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Figure A. 46 P (10 4(> C a 2 ), 1/4- Year Period Interior Motion. 




Figure A. 47 &P (10 4 ’C a 2 ), 1/4-Year Period Interior Motion. 
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Figure A. 48 Peak Enhanced and Reduced Cooling Values, Year 1- 




Peak Enhanced and Reduced Cooling Values, Year 2. 
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Figure A. 49 
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